Strong-Coupling Superconductivity in the Cuprate Oxide 
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Superconductivity in the cuprate oxide is studied by Kondo-lattice theory based on the t-J 
model with the electron-phonon interaction arising from the modulation of the superexchange 
interaction by phonons. The self-energy of electrons is decomposed into the single-site and multisite 
self-energies. It is proved by using the mapping of the single-site self-energy in the t-J model 
to its corresponding one in the Anderson model that the single-site self-energy is simply that of 
a conventional Fermi liquid, even if a superconducting order parameter appears or the multisite 
self-energy is anomalous. The electron liquid characterized by the single-site self-energy is a con- 
ventional Fermi liquid. The Fermi liquid is further stabilized by the resonating-valence-bond (RVB) 
mechanism. The stabilized Fermi liquid is a relevant unperturbed state that can be used to study 
superconductivity and anomalous Fermi-liquid behaviors. The so-called spin-fluctuation-mediated 
exchange interaction, which includes the superexchange interaction as a part, is the attractive 
interaction that binds d^2_j,2-wave Cooper pairs. An analysis of the spin susceptibility implies 
that, because of the electron-phonon interaction, the imaginary part of the exchange interaction 
has a sharp peak or dip at ±a;*, where lu* ~ ctJph in the normal state and ^ec ^ 5^0+ Wph 

in the superconducting state, where ojph is the energy of relevant phonons and ec is the su- 
perconducting gap. If the imaginary part has a sharp peak or dip at ±tj*, then the dispersion 
relation of quasi-particles has kink structures near ±cij* above and below the chemical potential, 
the density of states has dip-and-hump structures near ±aj* outside the coherence peaks in the su- 
perconducting state, and the anisotropy of the gap deviates from the simple ci^2_j^2-wave anisotropy. 

Keywords: strong coupling, high Tc, superconductivity, superexchange interaction, electron-phonon interac- 
tion, cuprate oxide, Fermi liquid, RVB, Kondo lattice 
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1. Introduction 

Although many experimental and theoretical studies 
have been performed since the discovery of the cuprate- 
oxide superconductor in the mechanism 

for high-temperature (high-Tc) superconductivity in the 
cuprate oxide is still contentious. Recently, Fe-based 
superconductors have been discoveredi^ Since both 
cuprate-oxide and Fe-based superconductors are in the 
vicinity of the Mott metal-insulator transition, it is cer- 
tain that strong electron correlations play a crucial role in 
the mechanism for high- Tc superconductivity. The most 
crucial issue is whether the normal state is actually an 
exotic Fermi liquid such as the resonating-valence-bond 
(RVB) stated* or it is nevertheless a conventional Fermi 
liquid. 

One of the simplest effective Hamiltonians for an elec- 
tron liquid in the vicinity of the Mott transition is the 
Hubbard model. According to Hubbard's theoryfi^iii 
when the on-site repulsion U is so large that U > W, 
where W is the bandwidth, the Hubbard gap opens be- 
tween the upper and lower Hubbard bands. According to 
Gutzwiller's theory ,^^'^'^i"'^'^ together with the Fermi-liquid 
theory^ii^ a narrow quasi-particle band appears near 
the chemical potential; this band is called the Gutzwiller 
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band in this paper. One may speculate that the density 
of states has a three-peak structure, with the Gutzwiller 
band between the upper and lower Hubbard bands. The 
approximations used in Hubbard's and Gutzwiller's the- 
ories are called Hubbard and Gutzwiller approximations, 
respectively, both of which are within the single-site ap- 
proximation (SSA). According to another SSA theoryjil 
the Gutzwiller band appears at the top of the lower Hub- 
bard band for n < 1, where n is the electron density per 
unit cell, which implies that it appears at the bottom of 
the upper Hubbard band for n > 1. The SSA that consid- 
ers all the single-site terms is rigorous for d -\-oo within 
the restricted Hilbert subspace where no order parameter 
exists fi^ii^'^i^ where d is the spatial dimensionality; this 
SSA is called the supreme single-site approximation (S"^ A) 
in this paper. The S^A theory is reduced or mapped to a 
problem of self-consistently determining and solving the 
Anderson model j^^'^'^i^^'^^ which is an effective Hamilto- 
nian for studying the Kondo effect. The three-peak struc- 
ture corresponds to that in the Anderson model, with the 
Kondo peak between two subpeaks. The Kondo effect 
has relevance to electron correlations in the vicinity of 
the Mott transition. The S'^A is also formulated as the 
dynamical mean-field theory^SiSe^^ (DMFT) and the 
dynamical coherent potential approximation^^ (DCPA). 
These three formulations are exactly equivalent to each 
other. 

According to Brinkman and Rice's theory,^" which also 
uses the Gutzwiller approximation, the Gutzwiller band 
vanishes when n — 1 and U > Ubr, with Ubr — W. Clus- 
ter DMFT (CDMFT)^i22i33^ j^a,s been proposed as an 



2 



extension of DMFT. According to the numerical results 
of BMFT^iSl'M and CDMFT^i^i^s^ ^^^^^ n ~ 1 and 
U > Ubr the Gutzwiller band appears to vanish and the 
ground state within the restricted Hilbert subspace ap- 
pears to be a Mott insulator. The numerical DMFT and 
CDMFT are consistent with Brinkman and Rice's theory, 
which treats the case of n = 1 , but are inconsistent with 
Gutzwiller's theory for n ^ I, which predicts that the 
ground state is a metal for n ^ I. 

In a previous paper,— it has been proved that the 
ground state in the S^A is a conventional Fermi liquid 
except when n — 1 and U/W — +oo, at least if an elec- 
tron reservoir is explicitly considered in the grand canon- 
ical ensemble. It is surprising that the ground states arc 
different between DMFT and the S'^A, which are equiv- 
alent to each other. This discrepancy is presumably be- 
cause the electron reservoir is not explicitly considered 
in the numerical DMFT, while it is explicitly considered 
in the proof. '^^ According to the proof, the conventional 
Fermi liquid constructed in the S'^A is a relevant unper- 
turbed state that can be used to study anomalous Fermi- 
liquid behaviors and ordered states in the vicinity of the 
Mott transition. A perturbative theory starting from the 
unperturbed state is simply Kondo-lattice theory i^^i^^'^"^ 
in which local electron correlations are, in principle, ex- 
actly treated in the S'^A, intersite electron correlations are 
treated by a conventional perturbation, and the sponta- 
neous appearance of an order parameter is treated by an 
anomalous perturbation. 

Although a multiband model is needed to explain the 
overall and precise features of both the cuprate-oxide and 
Fe-based superconductors, a single-band model can be ap- 
proximately used, at least for the cuprate oxide. For ex- 
ample, an effective Hamiltonian for the cuprate oxide is 
the d-p model, which considers d orbits on Cu ions and 
p orbits on O ions on Cu02 planes. When the on-site 
U is sufficiently strong, the d-p model is approximately 
mapped to the t-J model?^ which is a single-band model. 
The t-J model is also derived from the Hubbard model. 
The t-J model is one of the simplest effective Hamilto- 
nians for an electron liquid in the vicinity of the Mott 
transition or in the cuprate oxide. 

A mechanism for high- Tc superconductivity in the 
cuprate oxide was proposed in 1987>22ithe condensation of 
d2;2_j,2-wave Cooper pairs between heavy quasi-particles 
by the superexchange interaction. Since the supcrex- 
change interaction is as strong as J = —(0.10-0.15) eV 
between nearest neighbors;^ it is certain to play a role in 
the binding of d2.2_y2-wave Cooper pairs in the cuprate 
oxide. On the other hand, there is evidence that the 
electron-phonon interaction is strong in the cuprate ox- 
ide: the softening of the half-breathing modes near 
(±l,0)(7r/a) and (0,±l)(7r/a) in the two-dimensional 
Brillouin zone;^'^'^'^'^ where a is the lattice constant 
of the Cu02 planes, the softening of Cu-0 bond stretching 
modes near (±1/2, 0) (tt/o) and (0, ±1/2) (7r/a)^iil and 
kinks in the dispersion relation of quasi-particles42i^i^ 
A dip-and-hump structure outside the coherence peaks in 



the density of states has also been observed in the super- 
conducting state by tunnelling spectroscopy;^ A similar 
structure is observed in a conventional Bardeen-Cooper- 
Schrieffer (BCS) superconductor and is regarded as ev- 
idence that the binding of Cooper pairs is mainly due 
to the conventional electron-phonon interactionj ^^i^^i^'* 
However, the conventional electron-phonon interaction 
arises from the charge-channel interaction; thus, it can 
never be strong in the vicinity of the Mott transition 
because charge fluctuations are suppressed there. If the 
electron-phonon interaction is strong in the cuprate ox- 
ide, it must arise from the spin-channel interaction such 
as that arising from the modulation of the superexchange 
interaction by phononsi^i^ It will be interesting to study 
the role that the strong electron-phonon interaction plays 
in the mechanism for high- Tc superconductivity. 

One of the purposes of this paper is to formulate a 
theory of superconductivity based on the t-J model with 
the electron-phonon interaction. The other purpose is to 
apply it to superconductivity in the cuprate oxide. The 
preliminaries are given in §[2l The Kondo-lattice theory 
of superconductivity is formulated in §[31 In §[51 the for- 
mulation is applied to strong-coupling superconductivity 
in the cuprate oxide. A discussion is given in §[51 and a 
conclusion is given in §[6l An inequality that is crucial 
in the formulation in §[21 is proved in Appendix [XI The 
dynamical spin susceptibility of the t-J model is stud- 
ied in Appendix [B] The possibility of the spontaneous 
appearance of antiferromagnetic moments in the super- 
conducting state is examined in Appendix [Cl 



Preliminaries 



2.1. 



Electron-phonon interaction 



In the d-p model, according to field theory, the superex- 
change interaction arises from the virtual exchange of a 
pair excitation of d electrons between the upper and lower 
Hubbard bands. When nonzero bandwidths of the Hub- 
bard bands are ignored, the superexchange interaction 
constant is given by 



J 



uy 



1 



U 



(2.1) 



between d electrons on nearest-neighbor Cu ions, where 
V is the hybridization matrix between nearest neighbor p 
and d orbits, and are the depths of d and p orbits, 
respectively, and U is the on-site repulsion between d elec- 
trons. Equation (|2.ip is exactly the same as that derived 
by the fourth-order perturbation in V. When nonzero 
bandwidths of the Hubbard bands are considered, J is 
smaller than the value given by eq. (|2.ip . 

We assume that the Hamiltonian for phonons in the 
cuprate oxide is given by Hph = J2xq^)^q{b\gbxq + 1/2), 

where and b\q are the creation and annihilation op- 
erators, respectively, of a phonon with mode A and wave- 
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number vector and toxq is the energy of the phonon. 
The displacements of the ith Cu ion and the [ij]th O ion, 
which hes between the nearest-neighbor ith and jth Cu 
ions, are given by 



and 



hv. 



(2.2) 



, (2.3) 



where Vd.Xq and Vd^xq are real, Ri and -R[ij] = {Ri+Rj)/2 
are the positions of the ith Cu and [ij]th O ions, Md and 
Mp are the masses of Cu and O ions, respectively, N is 
the number of unit cells, and e\q — (eAq,^, exq.yi e\q,z) is 
the polarization vector. 

It is convenient to define a dual-spin operator by 

-Prig) = i^77ir(g') [S{q' + ^q)-S{-q'+^q)] , (2.4) 



q' 



where 



ii3)cr' 



(2.5) 



Here, cr = {(Tx, <Jy, <Jz) is the Pauli matrix, and dfco- 
are the creation and annihilation operators of d electrons, 
respectively, and ?7ir('7) are form factors defined by 



and 



Visiq) = cos(ga;a) + cos{qya), 



Vid{q) = cos((7^a) - coa{qya). 



(2.6) 



(2.7) 



It is assumed that the x and y axes are within Cu02 
planes and that the z axis is perpendicular to the Cu02 
planes. Two types of electron-phonon interaction arise 
from the modulation of J by the vibrations of O and Cu 
ionsj^i^ which are denoted by Hp and Hd, respectively: 



and 

'7Ad(9) = 2 eAq,:Esin (^^^ - CAq.ySin (^^^ . (2.10) 

Since the antiferromagnetic J is present between near- 
est neighbors, spin fluctuations develop near Qm = 
(±1, ±l)(7r/a).'^^ We denote the center wave number of 
spin fluctuations by Qsi — Qm + ^QsI, where |AQsf |a ^ 
1. According to eqs. (|2.4p and (|2.8p . two modes of 
spin fluctuations with wave number Qsf couple with a 
mode of phonons with wave number Qph — '^Qsf ± G 
or Qph = 2A(5st, where G is a reciprocal lattice vector. 
When spin fluctuations develop near Qm + AQsf, there- 
fore, phonons can be soft near 2A(5st.^- The electron- 
phonon interaction vanishes in the limit of IQphja ^ 
or jAQafja — > but it is nonzero for nonzero Qph or 
nonzero AQsf. When the electron-phonon interaction is 
strong, Qsi cannot be exactly Qm although Qsi — Qm, 
AQsf ^ but AQsf ~ 0. 

When coupled modes between spin fluctuations and 
phonons are sharp, 

Vsiq) = Vs{q ± G) = 5^ S (q' + iq) • S {-q' + \q) , 

(2.11) 

and Vdiq) — can be used as an approximation instead 
of eq. because ryis(O) = 2 and ?/irf(0) = 0. The total 
electron-phonon interaction, 7ioi-ph — Hp + TLd, is given 

by 



Tic 



Aq 



where Vs{q) is given by eq. I|2.1ip and 



Kx{q) = 



\j2LO\q 



Vxs{q) 



Cn 



'p,Xq 



/Mr, 



Vis 



CdVd, 



Aq 



(2.13) 



and 



'Hd 



'p,Xq 



5A-q 



X VXsiq) E '^ir {}q)Vriq), 

r=s,d 



hv, 



d.Xq 



^2NMdL0Xq 

X E Vxv{q)VT{q), 



5A-q 



(2.8a) 



(2.8b) 



where Cp and Cd are real constants, which are given in a 
previous paper 
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VXsiq) 



exq,x sm 



in(^)+eAq,,sin(^)], (2.9) 



2.2. Total effective Hamiltonian 
The t-J or t-J-Uoo model is defined by 



i(7 i^^jfy 



(2.14) 



where ta is the band center, Uia — dl^dia, Uj are transfer 
integrals, the summation (ij) runs over pairs of nearest- 
neighbor unit cells, and Si — {^/2)^^^i cr'^'' d\^dia' ■ 
The dispersion relation of bare electrons is given by 
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where Ri is the position of the zth unit cell. The band- 
width of E{k) is denoted by W. The on-site repulsion Uoc 
must be infinitely large to exclude any double occupancy. 
For the sake of convenience, Uoo is first treated as being 
finite and then the limit of Uoo/W — + +00 is taken in 
final results. It is assumed in this paper that the electron 
density per unit cell, n = I N)^^^ {riia) , is less than 
that at half filling such that < n < li^ 

Following a previous paper we explicitly consider an 
electron reservoir in the grand canonical ensemble. The 

reservoir is defined by Tires — J^ija^'ij^la-^j'^' ^^'^ 
finitesimally small but random hybridization between the 
reservoir and the t-J-Uoo model is defined by 



(2.16) 



where Aq = ±0^ is a nonzero but infinitesimally small 
numerical constant. In eq. (|2.16p . the summation (ij) 
runs over pairs of sites, the ith site in the t-J-Uoo model 
and the jth site in the reservoir, in the set TZ. Here, it is 
assumed that (( V[ij) )) = (( v*--^ )) = ^ind 



(2.17) 



where ((••■)) denotes the ensemble average for TZ and rih 
is the density of hybridization sites per unit cell, referring 
to the lattice on which the t-J-Uoo model exists. 
The total Hamiltonian to be considered is given by 

n = Ht-j + Hci-pi, + Hpu + n,es + V-nUt, (2.18) 



where /x is the chemical potential and A/t = J^icr Aa'^io + 
Sicr ^ia^i"- ^'^^ Sake of simplicity, it is assumed that 
the inversion symmetry exists in the total system aver- 
aged over the ensemble. 

When Uoo = 0, J = 0, and the electron-phonon inter- 
action is absent, the single-particle Green function for d 
electrons in the t-J-Uoo model averaged over the ensemble 
is given by 



iSn + E{k) - r(i£„) ' 



(2.19) 



Here, r(ie„) is the self-energy due to the hybridization 
with the reservoir. Since Ao = ±0"*", the second-order 
perturbation is sufficiently accurate to treat r(i£„) as 



r(ie,: 



1 



+ ^l-Eb{ky 



(2.20) 



where is the number of unit cells in the reservoir and 
i?f,(fe) ^ [l/Nb) E^^. expjifc • [R!^ - H;.)], where Rr is 
the position of the ith unit cell in the reservoir. It is 
assumed that no gap opens in the reservoir at the chemical 
potential: 



The electron number within the t-J-Uoo model is never a 
constant of motion in the presence of the electron reser- 
voir; thus, in general, the quantum-mechanically averaged 
number of electrons is a non-integer or an irrational num- 
ber. The presence of the reservoir with iV^ — *■ -\-oo or 
eq. (|2.2ip ensures the quality of the grand canonical en- 
semble that the averaged number of electrons is a contin- 
uous function of the chemical potential. 



2.3. Fermi-surface condition 



The Anderson model is defined by 



ka <7 

+ 7^ E (^'^4.'^- + VkdW.) , (2.22) 



where — d\da and Nc is the number of unit cells. 
The Green function for d electrons is given by 

I / / \ A I ^' I ^ 



lEn + p. - ed- ^^(ien) - - jde'^^^^ ^ 



(2.23) 

where Jl is the chemical potential, Y,„{ien) is the self- 
energy, and 

= ]^ E l^-^l'^e + M - E^ik)] . (2.24) 

k 

The Fermi surface of conduction electrons is defined by 
Ec{k) = fl. It exists when 

Aa(0) > 0, (2.25) 

is satisfied. This condition is called the Fermi-surface 
condition in this paper. 

The s-d model is another effective Hamiltonian for 
studying the Kondo effect. According to Yosida's per- 
turbation theory^*^ and Wilson's renormalization-group 
thcory,^^ the ground state is a singlet or a conventional 
Fermi liquid except when Js-d = 0, where Js-d is the 
s-d exchange interaction. Since the s-d model is derived 
from the Anderson model, the result for the s-d model 
implies that the ground state of the Anderson model 
is also a conventional Fermi liquid. The Bethe-ansatz 
solution for the Anderson model confirms that the 
ground state is a conventional Fermi liquid, at least 
when AA(e) is constant fS^iS'^lies jj-^ general, the nature 
of the ground state depends only on relevant low-energy 
properties, such as Aa(0), and high-energy properties 
only quantitatively renormalize the ground state, as 
demonstrated by renormalization-group theories for 
the s-d modelj^i^ When the Fermi-surface condition 
(|2.25p is satisfied, therefore, the ground state of the 
Anderson model is a conventional Fermi liquid except 
for the case corresponding to Js-d = 0, i.e., except when 
f7oo/AA(0) = -1-00 and Ud = {nd-\ + ridi) = 1. 



Imr(+iO) = 0+ > 0. 



(2.21) 
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Formulation 



3.1. Fermi liquid as an unperturbed state 



paper. When Uoo is nonzero but finite, J 7^ 0, or the 
electron-phonon interaction is present, the single-particle 
Green function for electrons is given by 



Since the superexchange interaction J is antiferromag- 
netic, only singlet superconductivity is studied in this 



J 



iSn - E{k) + /i - T.„{ien,k) -Ao-(ie„, k) 

-A;(ie„, k) isn + E{-k) - ^ + S_<,(-ie„, -fc) 



r 



(3.1) 



in the Nambu representation, where = 2TrkBT (n + ^) 
is the fermionic energy, T,cr{i-£n,k) is the self-energy of 
electrons, which includes r(i£„) defined by eq. (|2.20p . and 
Aa-{ien,k) is the superconducting order parameter; the 
Green function for phonons is given by 



[(iuji)'^ - ujI^]/ (2ujxq) - SA(iw£, q) ' 



(3.2) 



where uji = 2'!T£kBT is the bosonic energy and "Exlicui, q) 
is the self-energy of phonons. Because of the inversion 
symmetry, E{k) = E{-k), T,^{ien,k) = So.(ie„, -fc), 
Ao-(ie„, k) = Acr(i£„, -fe), uj\q = uj\-q, and T,x{iuji, q) = 
Yixiiuji, ~q). The determinant of Q^^{ien, k) is given by 



S^(ie„, k) = [i£„ - E{k) + /i - E^(i£„, fc)] 

X [ie„ + E{k) - ^1 + E_^(-i£„, A;)] 
- |A^(i£„,fc)|^ 



(3.3) 



Then, the Green function (|3.1|) is also described in such 
a way that 



t/o-(i£n,fc) 
where 

Go-(i£„, k) - 
and 



G'^(i£„,A;) iv(i£n,fe) 
F;(i£„,fc) -G-CT(-i£„,fc) 



i£„ -t- E{-k) - /X + S^(-i£„, -fc) 



A<,(l£„,fc) 

F^{ien,k) = — ^ — . 



(3.4) 



(3.5) 



(3.6) 



Feynman diagram in the site representation are classi- 
fied into single-site and multisite ones. Vertex corrections 
due to the random hybridization V can be ignored because 
they are 0(Aq), with Aq = ±0"*". The Green function for 
electrons in the site representation is given by 



G.,.(i£«) = ^5]e''=-(^'-^^-^G,(ien,fc), 



(3.7) 



for the diagonal component and 

F.jA^e^) - ^ ^e"=-(«-«^)F,(i£„,fc), (3.8) 



for the off-diagonal component. If only the lines of the 
on-site Uoo and the site-diagonal Giicr{i£n) appear in a 
diagram, the diagram is a single-site one. If a line of 
the intersite J, the site-off-diagonal Gijcr{i£n) with i ^ j, 
the phonon Green function, or the off-diagonal Fjjo-(i£n) 
appears in a diagram, the diagram is a multisite one. Ac- 
cording to this classification, the self-energy of electrons is 
decomposed into the single-site So-(i£n) and the multisite 
AEo-(i£n, k) in such a way that 



E^(i£„, fc) = S^(i£„) + AS^(i£„, fc) + r(i£,; 



(3.9) 



In the Anderson model, all the Feynman diagrams are 
simply single-site diagrams. It is assumed that — /i and 
Uoo of the Anderson model are exactly the same as — 
and Uoo of the t-J-U 00 model, respectively: ed — fi — Cd — M 
and Uoo = Uoo- The same Uoo appears in the Feynman 
diagrams of the Hubbard and Anderson models. If Aa(£) 
is determined to satisfy 



Giia{^£n) — G(j(i£„), 



(3.10) 



the single-site self-energy is equal to the self-energy for 
the Anderson model: 'Sg-i^Sn) = So.(i£„). It follows from 
eq. ([XTU)) that 



AA(e)=Im S,(e + iO) + G,^i(e + iO) 



(3.11) 



Equation ([3T0)) or ((3ll|) is the mapping condition to the 
Anderson model. In general, AA(e) depends on T, i.e., 
the mapped Anderson model itself depends on T. 

According to eq. p.lOp . the electron density and the 
density of states of the t-J-Uoo model are the same as 
those of the Anderson model: n — Ud and 

p{e) = --l■mG^^„{e + iO) = -ilmG^(e + iO). (3.12) 



The mapping condition p. lip is iteratively treated 
to obtain the eventual self-consistent AA(e); not only 
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AA(e) but also ^^(e), AEo.(e + iO, fc), Ao-(ie„,fc), and 
S>(aj + iO, k) should be self-consistently calculated to sat- 
isfy eq. (|3.1ip . It is proved in Appendix [A] that 



AA(e) > -Imr(e + iO). 



(3.13) 



According to eqs. (|23T|) and ((3?T3)) . Aa(0) > 0, i.e., the 
Fermi surface condition ()2.25|) is satisfied at each step of 
the iterative process. Unless Uoo/W = +oo and n — 1 
exactly, therefore, the ground state of the Anderson model 
is a conventional Fermi liquid, so that the single-site self- 
energy for the ground state is that of the conventional 
Fermi liquid even if the order parameter is nonzero or the 
multisite self-energy is anomalous. 

We consider the Anderson model in the presence of an 
infinitesimally small chemical potential shift A/iA and an 
infinitesimally small Zeeman energy h. The self-energy 
Scr(e + iO) for the Anderson model is expanded in such a 
way that 

S,(e + iO) = So + (l-0o)e 

+ 0(eVfcBTK) +0(fcBrVTK), (3.14) 

where So, 0c > 0, (j>c > 0, > 0, and Tk are all real. 
Here, Tk is the Kondo temperature and fceTk is the en- 
ergy scale of local quantum spin fluctuations for not only 
the Anderson model but also the t-J-Uoo model. In the 
following part of this paper, T <^ Tk or 



Tc « Tk, 



(3.15) 



is assumed, where Tc is the superconducting critical tem- 
perature and (j)s/4'c is denoted by Wg'- 



Ws = (f>s/(t>c 



(3.16) 



which is the so-called Wilson ratio. In general, 20c = 
<j)s -\- 4>c- When n ~ 1 and Uoo/W 1, local spin fluc- 
tuations are developed but local charge fluctuations are 
suppressed, so that (f)s ^ 1 and 0c <C 1. Then, 0e ^ 1 
and Ws ~ 2. 

When the order parameter Ao-(ien,fe) is ignored, the 
diagonal part of the Green function is given by 



Go-(ie„, k) 



1 



1 



>c ie« - Co(fc) - AE<^(i£„, fe) 
[incoherent term], (3-17) 



where 



and 



eo(fc) = ^[So + i?(fe)-/^], 

0c 



AS^(ie„,fc) = J-AS<,(i£„,fc) 



(3.19) 



The first term in eq. (|3.17p is the coherent term, which 
describes the quasi-particle band with bandwidth W* = 



0(fcB'7k) at the top of the lower Hubbard band. The 
second term in eq. p.l7p is the incoherent term, which 
describes the lower Hubbard band; the upper Hubbard 
band lies at an infinitely high energy in the limit of 
Uoo/W ^ oo. 

It is convenient to define other renormalized quantities 



A,(i£„, fe) = ^A^(ie„, fe), (3.20) 

0e 

S<j(ie„, fe) = [ie„ - ^o(fe) ^ Ai;^(ie„, fc)] 

X [ie„ + eo(fe) + AI]_^(-ie„, fe)] 
- |A,(i£„,fe)|', (3.21) 



G.(ie„,fe) = i^» + eoW + Ays„,fe) ^ ^3^22^ 



and 

When the incoherent part is ignored, the single-particle 
Green function is given by 



1 



C/c,(ie„,fe) = -=-t/cr(ien,fe) 



(3.24) 



where 



3.2. Intersite exchange interactions 

The polarization function in spin channels is also 
decomposed into the single-site TTs{iuJe) and multisite 



(3.26) 



The spin susceptibilities of the Anderson and t-J-U^ 
models are given, respectively, by 



(3.18) and 



1 - ?7oo7rs(iw^) ' 



1 - [jJsiq) + Uool 7rs(iw^,q) ' 



where 



Js{q) = j^'^J exp [iq ■ {R, - Rj)] 



(3.27) 



(3.28) 



(3.29) 
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FIG. 1: Feynman diagram corresponding to Js{-ph{ii^i,q)- A 
wavy solid line denotes Uoo- A hatched circle denotes the 
irreducible single-site vertex function in spin channels, and an 
arrowed solid line denotes the Green function for electrons; 
a unit composed of a hatched circle and two arrowed solid 
lines is the irreducible polarization function in spin channels. 
A hatched triangle denotes the electron-phonon interaction, a 
dashed line denotes the spin susceptibility, and a dotted line 
denotes the Green function for phonons. The i'th and j'th 
sites are nearest-neighbor Cu sites of the ith and jth Cu sites, 
respectively. The [ii'Jth and [jj'lth sites are O sites between 
the ith and i'th Cu sites and between the jth and j'th Cu 
sites, respectively. 



It should be noted that in the limit of Uoo/W +00, 

L/oo^^s(iw£) ^ 1, (3.30) 

and C/ooA7rs(iwf,q) O [l/UocXs{iuJi,q)]. A physical 
picture of Kondo lattices is that local spin fluctuations 
on different unit cells interact with each other by an in- 
tersite exchange interaction. Then, the exchange interac- 
tion, which is denoted by Is{iujf , q), is defined by 



1 - \ls{\uJi,,q)xs{T'^^i) 



It follows from this definition that 

Isiiuii,q) = Jsiq) + 2U^ATTsiiuji,q), 



(3.31) 



(3.32) 



in the limit of Uoo/W +00. 

The term 2U^ATTs{iuj£,q) in eq. (|3.32p is composed of 
several terms. One is an exchange interaction arising 
from the virtual exchange of a pair excitation of quasi- 
particles. When the reducible and irreducible three-point 
single-site vertex functions in spin channels are denoted 
by As(ie„-f iwf, ie„; itue) and As(ie„ ie„; iw^), respec- 
tively, it follows that 



As(ie„ + ia;£,ie„;iwf) 



As (ie,: 



1 - UooTfs{iuJi) 



(3.33) 



according to the Ward relation.^'' Since (j)s — As(0,0;0), 
it follows that 



(3.34) 



in the limit of £„ ^ 0, — > 0, and Uoo/W — > +00. When 
eq. (|3.34p is used as an approximation for small e„ and 



to I such as |e„| <C fceTk and \u!i\ <^ k^T^, the exchange 
interaction is given by 



2 r 



jQ{'^^e,q] 
where 
P{iuJi, q) = - 



P{iuje,q) - P{iuji) 



(3.35) 



' N ^ 

e„fe 



[a(i 



1£„ 



iw£, fc + q)Gaiien,k) 
F^{ien + iuje,k + q)F^{ien,k)] , (3.36) 



and 
-P(iw£) 



X! ^'T(i£« + iw£,fe)G'^(i£„,p), (3.37) 



iV2 



which is derived in the random-phase approximation 
(RPA) for a pair excitation of quasi-particles. In 
eq. (|3.35p . the single-site term P{\loi) is subtracted to ex- 
clude any double counting. Another term is an exchange 
interaction arising from the virtual exchange of a cou- 
pled excitation of spin fiuctuations and phonons, whose 
exchange process is schematically shown in Fig. [TJ When 
eq. (j3.30p is made use of, the exchange interaction is sim- 
ply given by 

Jsf_ph(it^^ q) = 2—— ^ Kl{q)Dx{iuJi + icj„, 2q ± G + p) 



N 

UJnP 

X Xs(iWn,9+P), 



(3.38) 



where K\{q) is given by eq. ()2.13p . When only these 
terms are considered, it follows that 

Is{iu}e,q) = Js{q) + JQ{iuJi,q) + Jsf_ph(iw£, q). (3.39) 



When eq. (|3.34p is used as an approximation, the spin- 
fiuctuation-mediated exchange interaction is given by 

J(C/A,)'[Xs(i^^q) - Xsi/y^i)] = i^^/:(ic^f,g), (3.40) 

where denotes As(0, 0; 0) and 

is{iujt,q) 



Isi^^e,q) = 



^(iw^,9)Xs(iwf) 



(3.41) 



In eq. (|3.40p . the single-site term is subtracted to ex- 
clude any double counting, and two 0^ appear as effective 
three-point vertex functions. According to eqs. (I3.40p and 
(|3.4ip , the spin- fluctuation-mediated interaction is simply 
the exchange interaction I*{iuje,q); the bare Is{i(-ue,q) is 
enhanced to I*{iLL>i, q) by spin fluctuations. 



3.3. Gap equation 

The enhanced exchange interaction I*(iu;e,q) includes 
no single-site effect from its definition; thus, it does not in- 
clude the strong effective on-site repulsion between quasi- 
particles that arises from single-site correlations due to 
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the infinitely large on-site Uoo- When the effective on-site 
repulsion is denoted by r^, the gap equation is given by 



N 

Ag(i££,p) 
S<^(ie£,p) 



k) 



(3.42) 



where ^^(ien, k) is defined by eq. (|3.2ip . Here, the energy 
dependence of is ignored because its energy scale is 
/cbTk, which is assumed to be much larger than k^Tc. 

When only the first-order term in I*{\LjJi,q) is consid- 
ered, for example, the multisitc self-energy of electrons is 
given by 



AS^(ie„,fe) 



N 



X /;(ief-ie„,fc-p)G„(ie<,,p), (3.43) 



where GaiiSmk) is defined by eq. (13. 22^ : T-^i is not in- 
cluded in eq. (|3.43p to exclude any double counting. If 
superconducting fluctuations develop, their renornializa- 
tion, which is of higher order in q), should be con- 

sidered. The self-energy of phonons T,x{iuje,q) can also 
be perturbatively calculated in 7ici-ph, as examined in a 
previous paper 

The theory in this paper is similar to the conventional 
theory of strong-coupling superconductivity except that 
the unperturbed state, which is a conventional Fermi liq- 
uid, should be self-consistently constructed with other 
quantities such as Ao-(ie„,fc), ATia-0-Smk), Y,\{iuji,q), 
I*{iuji,q), and so forth, to satisfy the mapping condi- 
tion (jsiU)) or ((3ll|) . 



4. Application to the Cuprate Oxide 

4.1. RVB stabilization mechanism 

The exchange interaction is expanded in the Fourier 
series in such a way that 

/;(iw,, q) = I^iiuje) + 2^ /;(ic^,)7,„,(q), (4.1) 



where 



Here, Rq and Rj are the lattice vectors of the origin and 
its jth nearest neighbor, respectively, and r]ns (q) is the s- 
wave form factor of the nth nearest neighbors: ?7os('Z) = 1, 
r]is{q) is defined by eq. (HJ]), ri2s{q) = cos[(g2: + qy)a] + 
cos[{qx — <ly)o]^ and so forth. In general, 

J+[/r(0)- J]ai(ic^,), j = l 



I*{iu;i) = i ' '\iY' .J"^v'--^' ^ ; (4.3) 



where aj{iuji) is an analytical function that satisfies 
aj{0) = 1 and [aj(iwf)] |^^|^_|_^ = 0. When the con- 
stant term J or the superexchange interaction Js{q) is 
only considered in eq. ()3.43p , the self-energy does not de- 
pend on energy; thus, it is simply denoted by AI]j(fc) 
here. The diagonal part of the Green function is given by 
Ga{ien,k) = 1/ [ie„ - £,{k)], where 

ak)^^[f:a + E{k)-fi]+Atj{k), (4.4) 



is the dispersion relation of renormalized quasi-particles 
to be determined. When T = K, Eq — M is determined 
so that it satisfy the Fermi-surface sum rule: 



(4.5) 



where H{x) is the Heaviside function defined by H{x > 
0) = 1 and H{x < 0) = 0. It follows that 

Atjik) = ^VK2jr7i,(fc)l^i/[-^(p)/W^]r;i,(p), 

p 

(4.6) 

at T = K. Since J is antiferromagnetic, the bandwidth 
of quasi-particles is broadened by this renormalization; 
even in the limit of (pc +oo, the bandwidth of quasi- 
particles is 0(|J|). The unperturbed Fermi liquid con- 
structed in the S'^A is further stabilized. The stabiliza- 
tion is due to the formation of an itinerant local singlet or 
a resonating valence bond (RVB) on each pair of nearest 
neighbors and is similar to that in the mean-field RVB 
theory,— although the stabilized liquid is simply a con- 
ventional Fermi liquid. The stabilized Fermi liquid is a 
more relevant unperturbed state than the Fermi liquid con- 
structed in the S'^A. 

The density of states for electrons at the chemical po- 
tential is given by 



p(0) = 



f:o + E{k)+4>,Atj{k)~ fi . (4.7) 



It is different from that for quasi-particles, which is de- 
fined by 



(4.8) 



Note that p(0) = O[l/(0e| J|)] in the limit of ^ +oo. 
The density of states p(0) is greatly reduced by the ex- 
istence of 0eAEj(fc) with <^e ^ 1, although p*(0) is 
not reduced unless J is large. According to the Fermi- 
liquid relationji^iiS the specific-heat coefficient is given 
by 7c = (2/3)7r^fc|p*(0). The quasi-particle bandwidth, 
which is defined by W* = W/(f)c, is estimated to be W* ~ 
l/p*(0) ~ 0.3-0.4 eV from the observed specific-heat 
coefficient)^ which is, for example, about 10 mJ/niol-K^ 
for optimal-doped cuprate oxides, in which is the high- 
est as a function of doping. Since W = 3-4 eV according 
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to a band calculation, W/W* ~ 10, which implies that 
~ 10. On the other hand, the superexchange interac- 
tion is as strong as J = —(0.10-0.15) eV4S It is certain 
that the RVB stabilization mechanism is crucial in the 
cuprate oxide with n ~ 1. It will be interesting to exper- 
imentally study the reduction of p(0) for electrons, not 
p*(e) for quasi-particles, in the cuprate oxide with n ~ 1 
to obtain evidence that the RVB stabilization mechanism 
is actually crucial. 

If the Kondo temperature is defined by fceTk = 
l/2p*(0), then /cbTk ^ 0.15-0.2 eV for optimal-doped 
cuprate oxides and fcelk > 0.2 eV for over-doped and 
under-doped cuprate oxides. The condition (j3.15p . Tc <^ 
Tk, is satisfied in any cuprate-oxide superconductor. 



4.2. Kmks m the normal state 

The exchange interaction /*(iw£, q) or its Fourier com- 
ponent I*{iuji) should be self-consistently calculated us- 
ing the self-energies of electrons and phonons, the order 
parameter, the spin susceptibility, and so forth. Since it 
is a difficult task to complete the self-consistent process, 
a phenomenological theory is developed in the following 
part of this paper; each crucial effect is studied in a non- 
self-consistent manner or independently by assuming that 
the analytical continuation of aj{iuj^) onto the real axis 
can be approximately given by 



where 



aj{ijj ± iO) 



2v, 



to -\- Vj ± i7j — jAj ± i7j 



(4.9) 



where Vj and 7j are positive constants. Figure UKa) shows 
aj{iji) -\- iO) for several ^jjvj. According to the analysis in 
Appendix [B] it is reasonable to assume that Vj ~ Wph 
in the normal state and that (l/2)eG ^ Vj ^ (l/2)eG + 
Wph in the superconducting state, where Wph is the energy 
of relevant phonons and ec — max[2|Acr(+iO, A;)|]. For 
the sake of simplicity, it is also assumed that the density 
of states for quasi-particles in the unperturbed state is 
constant, so that 



(4.10) 



Here, D* ~ W/2ipc is half the quasi-particlc bandwidth 
and p*o = 1/(21?*). 

Since the self-energy due to the on-site component does 
not depend on fc, it is simply denoted by ASo(ien) here; 
the retardation effect is crucial in the on-site effect. When 
Go.(ie„,fc) = l/[ie„-^(A;)] is used in eq. (|3.43|) . Ai;o(ie„) 
is given by 



AI]o(ie„) = .90 



J-D- 



aoiiei - ie„) 



l£i ■ 



(4.11) 



90 



(4.12) 



Since (1/iV) Eg Js{q) = and (1/iV) Eg Jq(}uj,, q) = 0, 
/q(0) is simply given by 



J,t-ph(0,g) + -/,2(0,q)x,(0,q) 



(4.13) 



It is obvious that /o(0) is positive and that 50 > 0. It 
should be noted that ReAi;o(e + iO) and ImAEo(e + iO) 
are odd and even functions of e, respectively, within the 
phenomenological model used in this subsection. 

In this subsection, vq/D* — 1/4 is assumed; results for 
low-energy parts are not sensitive to the value of / D* 
provided that v^jD* is sufficiently small. Figure [2jb) 
shows AEo(e-hi0)/(5oJ^o) for7o/j^o = 0.1, 0.2, 0.4, and 0.8. 
The renormalized dispersion relation of quasi-particles, 
which is denoted by f*(fc), is approximately determined 

by 



r (fc) = C(fc) + Re {AEo r(fc) + iO]} 



(4.14) 



FigurelUc) shows C*(fc) as a function of ^(fe). Since ^(fc) 
is a smooth function of fc, ^*(A;) shows kinks as a function 
of fc at approximately iz^o above and below the chemi- 
cal potential. Figure [2jd) shows the density of states for 
renormalized quasi-particles: 



P5 



Im / de' -. 

TT Ld' e-e'- AEo(e + iO) 



(4.15) 



When a kink appears at approximately e ~ ±fo in the 
dispersion relation ^*(fc), a small structure appears at 
approximately e ~ -^vq in the density of states p*(e). 

4.3. Dip-and-hump structure in the superconducting state 



When the on-site and the nearest-neighbor /*(i(jjf) 
are only considered, only s and dx2-y2 waves are possi- 
ble. Since T^i is strongly repulsive and r]'lj^{k) <C 1 on 
the Fermi surface for n ~ 1, when n ~ 1, T,, of the s 



wave is much lower than T„ of the d^ 



Thus, 



the c?j.2_j,2 wave is only considered here. Since Ff^ has no 
effect on the d^2_y2 wave, the order parameter A^{ien, k) 
is decomposed into the form factor 7]id{k) and an energy- 
dependent part, which is called a gap function and is sim- 
ply denoted by AG(ien) here: 



A^(i£„,A;) = -77id(fc)AG(i£„). 



(4.16) 



To simplify the numerical processes, the renormalization 
of quasi-particles studied in ii l44.2l is ignored, and 



r?id(fc) = 2cos(2(/3), 



(4.17) 
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FIG. 2: (a) Phenomenological function aj{uj + iO), (b) self-energy AEo(£ + iO) for D* /vq — 4, (c) kink structure for go — 0.3 
and D* /uo — 4, and (d) density of states p*(e) for D* /vq = 4 and three values of go: (i) go = 0.3, (ii) 0.5, and (iii) 0.7. In each 
figure, (red) 'jj/uj = 0.8, (yellow) 0.4, (green) 0.2, and (cyan) 0.1. In (a) and (b), the solid and dashed lines show the real and 
imaginary parts, respectively. In (c), the relation ^*(fc) = f(fc) is shown by a black line for comparison. In (d), p*(e) for go = 
is shown by a black line for each value of go for comparison. 



where ip — taii^^ {ky / kx) , is used as an approximation 
instead of eq. ()2.7|) . The particle-hole symmetry does not 
exist in the original t-J model, in general. However, the 
particle-hole asymmetry is not crucial here, similarly to 
in the BCS theory. Thus, the constant and symmetric 
density of states given by eq. (|4.10p is assumed. When 
the summation along the imaginary axis is transformed 
into the integration along the real axis in eq. (I3.42p . the 
gap equation is given by 



AG(e + iO) 



-gi / d( 

J-oo Jo 
J-oo \ 



dip 



2-K 



2cos^(2^) 



Im 



/r(z + io) 



AG(2 + e + iO) 
S(z-He + iO,^,V5) 



iO) 



/r(o) 



Im 



AG(^ + iO) 
S(z + iO,^,(p) 



>, (4.18) 



at r = K, where 
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(4.19) 



is the dimensionless coupling constant and 



S(e-fi0,^,^) = (e + i0)2 



+cos2(2v3) |AG(e-l-iO) 

(4.20) 



In eq. (j4.18p . the integration over ^ is extended from 
—D* < C < +D* to — oo < ^ < +00 to simplify the nu- 
merical processes. Instead of this, a cutoff is introduced 
into the superexchange interaction such that the constant 
J is replaced by Jaj{uj -f iO), where 



aj{uj + i{)) = H{l-\uj\iej). 



(4.21) 



Here, Oj is the phenomenological cutoff parameter; 
vi/Oj = 1/A oi Oj/vi = 4 is assumed. Although the 
analyticity of the summand or integrand is assumed in 
the transformation from the summation along the imagi- 
nary axis to the integration along the real axis, this cutoff 
function is not analytical. To confirm the relevance of the 
cutoff model of eq. (|4.21|) . another cutoff model is also 
studied: aj{uj + iO) ~ Hjj{u;)/H.yj{0), where 



H^j{uj) = — In 



U! -9j + ijj 



(4.22) 



is an analytical function that satisfies [Hjj(uj)]^^_^^q = 
H{1 — \uj\/6j). Here, jj = O.Wj is assumed to re- 
move the logarithmic singularity in ImiJ^^ (w) . These 



11 



TABLE I: Adjusted gi for the analytical (A) and non- 
analytical (non-A) models with ji/vi = 0.1 to reproduce 
2Ag(0)/i/i = 1 or 2Ag(0)/6»j = 1/4. 



ri 


A 


non-A 


0.3 


0.5455 


0.6080 


0.5 


0.4745 


0.5175 


0.7 


0.4323 


0.4605 


1 


0.3893 


0.3928 



two models are called the non-analytical and analyti- 
cal models, respectively; the main difference between 
them is that Imaj(u; -I- iO) = in the non-analytical 
model but linaj{uj + iO) 7^ in the analytical model. 
When J/W < 0, l*{0,q) is antiferromagnetic or l*{0,q) 
has its maximum near Qm = (ilj il) (7r/a). Then, 
Il{0)/W < so that gi > 0. When a constant ri is 
defined by n = J//* (0)', + iO)//i*(0), which appears 
in the gap equation (|4.18p . is given by 

I^ioj + iO)//i*(0) = riaj{uj + iO) + (1 - ri)ai(w + iO). 

(4.23) 

Since the phase of AG(e + iO) is arbitrary, it is chosen 
in such a way that AG(+iO) is real and positive; thus, 
it is simply denoted by Ag(0) here. Since the particle- 
hole symmetry exists in the simplified phenomenological 
model, the real and imaginary parts of AG(e+iO) are even 
and odd functions, respectively. Figure [3] shows AG(e+iO) 
as a function of e for the two models with ji/vi = 0.1. 
To compare the non-analytical and analytical models, the 
dimensionless coupling constant gi is adjusted such that 

2Ag(0)/z/i = 1, (4.24) 

or 2Ag(O)/0j = 1/4. The adjusted values of gi for 
= 0.1 are shown in Table IH The gap function 
|AG(e + iO)| is larger in the analytical model than in the 
non-analytical model because of ImQ;j(a; ± iO), but no 
crucial difference can be seen between them in the low- 
energy region of |e| < (1.5-2.0)2Ag(0). Results for only 
the non-analytical model are shown in the following part 
of this paper. The adjustment of eq. (|4.24p for gi is made 
in each case. 

Figure [4] shows the density of states for Bogoliubov's 
quasi-particles, which is given by 

J_oo Jo 27r V Ti-y s(£-Hi7q,C,</') 

(4.25) 

where ^q/Oj — > +0 is assumed here. The coherence peak 
appears at e ~ Ag(0) for ri > 0.2. A dip and a hump 
appear at e ~ \yi + Ag(0)] /2 and e ~ i^i, respectively, 
for 0.2 < ri < 0.3. For ri = 0.19, on the other hand, three 
peaks appear, which implies that a first-order transition 
or a sharp crossover occurs between ri > Tc and ri < Tc, 



(iv) 




£/2Ag(0) 



FIG. 3: Gap function AG(e + iO) for Oj/ui = 4 and 
7i/z/i =0.1: (solid line) the non-analytical model and (dashed 
line) the analytical model; (green) Re [Ac (e -I- iO)] , (red) 
Im [AG(e + iO)], and (blue) |AG(e + iO)|; (i) n = 1, (ii) 
ri = 0.7, (iii) r\ = 0.5, and (iv) ri = 0.3. In (i), the imaginary 
part is vanishingly small in the non-analytical model. 



(V) 




e/2A^(0) 



FIG. 4: Density of states p*(e) for 9,j /vi — 4 and "/i/i'i = 0.1: 
(i) n = 0.7, (ii) n = 0.5, (iii) n = 0.3, (iv) n = 0.2, and 
(v) ri = 0.19; for comparison, p*(e) for ri — 0.7, 0.5, 0.3, and 
0.2 are also shown in (v). As can be seen in (v), p*(e) hardly 
depends on ri in the low-energy region of |e|/2AG(0) < 0.25 
or |e|/AG(0) < 0.5 because Ag{0)/i^i ~ 0.5 is assumed in 
every case, but the apparent size of the gap for r = 0.19 is 
about twice as large as that for < ri < 0.2. In the case 
of ri — 0.19, the sum rule for the density of states p*(£) is 
significantly violated because of non-analyticity involved in the 
non-analytical model; thus, the sharp three-peak structure for 
ri = 0.19 should not be seriously considered and should be 
regarded as an approximate or qualitative one. 



with Vc — 0.2. The transition or crossover is studied in 
the next subsection of ? 14 4.41 Figure [5] shows the density 
of states for ri — 0.2 and three values of ^ijvi: 0.1, 0.2, 
and 0.4. The dip-and-hump structure can only be seen 
for sufficiently small 71/2^1. 

According to the analysis in Appendix [Bl it is likely 
that Ag(0) < 1^1 < Ag(0) -l-Wph. The observed dip-and- 
hump structur^i is consistent with results in this paper. 
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(iii) 




e/2Ag(0) 



FIG. 5: Density of states p*(e) for 9j/vi = 4 and ri = 0.2: 
(i) "li/vi = 0.4, (ii) 0.2, and (iii) 0.1. 

4.4. Small-gap and large-gap phases 

The dispersion relation of quasi-particles in the super- 
conducting state, which is denoted by £,±{k), is also de- 
fined by the pole of the Green function. It is approxi- 
mately given by a solution of 

c±{k) = ±^je{k) + ^^uk)\Ko[e±{k) + -M\\ 

(4.26) 

According to this equation, an effective gap as a func- 
tion of If, which is denoted by eG(v), is also defined by a 
solution of 

\eG{v) = |cos(2(^)Ag [iecM + iO] | , (4.27) 

which is the solution of eq. (I4.26P for k satisfying ^{k) — 
0; eqs. (|4.10p and (|4.17p are assumed. As is implied by 
Fig. [6l there is only a single solution when ri < 0.2 while 
there are three solutions when ri = 0.19 and \if \ is small. 
Because of the appearance of extra poles in the Green 
function, the gap equation gives a larger gap function 
when ri < rc — 0.2 than when ri > rc — 0.2. The jump 
between ri > Tc and ri < rc is that between a small- 
gap phase for ri > rc and a large-gap phase for ri < rc- 
Within the numerical treatment of this paper, the jump 
appears to be a first-order transition rather than a sharp 
crossover. 

Since AG(e + iO) depends on e, the (^-dependence of the 
effective gap cq {f) defined by eq. (|4.27p deviates from the 
cos(2(^) dependence, as shown in Fig. [T] The deviation 
of the gap anisotropy from the simple c?-wave anisotropy 
has been actually observed in the cuprate oxide>^ 

Figure [8] shows the density of states when the lifetime 
width 7q of Bogoliubov's quasi-particles is considered in 
eq. (|4.25p . The two low-energy peaks, which appear for 
ri = 0.19 and 7q -1-0, are not resolved even for small 
7q/2Ag(0) = 0.03; the resolution is worse if the life-time 
effect is considered in the gap equation. It is likely that 
the two low-energy peaks are not resolved in the cuprate 
oxide even if the large-gap phase actually appears. 



(iii) 




e/2A^(0) 



FIG. 6; Gap function |AG(e4-iO)| for Oj/ui = 4 and 71 ///i = 
1: (i) ri = 0.3, (ii) 0.2, and (iii) 0.19. The relation |AG(e -I- 
iO)| = e is shown by a dotted line. 



o 




2ipln 



FIG. 7: Anisotropy of the effective gap eG('p) for 9,j/v\ — 4 
and 7i/i^i = 0.1: (dotted line) ri = 0.3, (dashed line) 0.2, and 
(dot-dashed line) 0.19. For comparison, cos(2(p) is shown by a 
solid line. Note that £0(^5) is multi-valued for ri = 0.19; it is 
expected in such a case that the observed gap increases much 
more rapidly as ip decreases than cos(2(p) does. 



(iii) 




12 3 

FIG. 8: Density of states p*(e) for 9j /u\ — 4 and 'yi/fi — 0.1: 
(i) n = 0.3, (ii) 0.2, and (iii) 019; (solid line) 7q/2Ag(0) = 0, 
(dotted line) 0.03, (dashed line) 0.1, and (dot-dashed line) 
0.3. The sum rule for p*(e) that is significantly violated in 
the case of ri — 0.19 is approximately recovered when the 
phenomenological lifetime width 7q is introduced. 
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5. Discussion 

The exchange interaction + iO,q) given by 

eq. (|3.4ip is composed of various terms: 

+ iO, q) = Is{uj + iO, q) + i/f (w + iO, q)xsii^ + iO, q), 

(5.1) 

where Is{uj + iO, q) is given by eq. p.39p or Is{uj + 
iO, q) = Js{q) + Jq{uj + iO g) + Jst-ph(^^ + iO, g) and 
Xs(w+iO, g) is given by eq. (|3.3ip orxs(w+i0, g) = 
iO)/ [1 - (l/4)/s(w + iO, g)xs(w + iO)]. A crucial issue is 
to determine the main attractive interaction among the 
various terms that binds dx^ -y^-^tcve Cooper pairs. Ac- 
cording to Appendix iBl Im[Jsf_ph(w + iO, g)] has a sharp 
peak at approximately Wph in the normal state and at 
approximately Ag(0) -f-Wph in the superconducting state; 
when Jsf_ph(ti^ + iO, g) has such a sharp peak, it is proba- 
ble that Xs(w -f iO, g) also has a similar structure in the 
same or a similar energy region. According to the study 
in 14 4.31 the appearance of the dip-and-hump structures 
outside the coherence peaks can be explained by the exis- 
tence of such a peak in the imaginary part of /* (w -I- iO, g) . 
If the dip-and-hump structure appears, 

Jsf-ph(w + iO, g) + i/2(cj + iO, q)xs[i^ + iO, g), (5.2) 

must be at least as effective as the superexchange in- 
teraction Js{q)', if no dip-and hump structure or only a 
small one appears, the superexchange interaction Js{q) 
must be the main attractive interaction. It is certain that 
Xs{i^ + iO, g) is enhanced in the vicinity of the Neel state. 
As studied in Appendix [Bl Jsf-ph{i^ + iO, g) includes the 
convolution of Xs('^ + iO, g) and the Green function for 
phonons. Therefore, it is reasonable that the dip-and- 
hump structure is larger in under-doped cuprates than in 
over-doped cuprates. 

The kink structure in the dispersion relation is mainly 
due to the on-site component of I*{lu + iO, g), and the 
dip-and-hump structure in the density of states is mainly 
due to the nearest-neighbor component of I*{u! + iO, g). 
Although their appearances are closely related with each 
other, they may be different from each other; one may 
be observed but the other may not be observed, or their 
characteristic energy scales may be slightly different from 
each other such as vq ^ vi and vq ~ vi. 

Three types of electron-phonon interaction are possible 
within the t-J model: one arising from the modulation 
of the band center or site energies by phonons, one from 
that of transfer energies ti^j by phonons, and the one 
considered in this paper. Since the first type of interac- 
tion couples with charge fluctuations, the single-site ver- 
tex correction divided by the mass enhancement factor is 
given by <^c/^e- Since 1 and ^ 1, the first type 

can never play any role in the t-J model with almost half 
filling of electrons. In the second type, no single-site ver- 
tex function can appear and l/4>c appears in every line 



of the electron Green functions. The second type cannot 
play a crucial role in the t-J model. Only the third type 
can play a crucial role in the t-J model with almost half 
filling of electrons. 

If the electron-phonon interaction is ignored, the cen- 
ter wave number of antiferromagnetic spin fluctuations 
should be exactly Qm = (±1,±1) (vr/a) or very close to 
Qm- When the electron-phonon interaction is strong, 
however, it cannot be Qm- The electron-phonon in- 
teraction can assist the development of antiferromag- 
netic spin fluctuations at wave numbers slightly differ- 
ent from Qm and an antiferromagnetic order with such 
wave numbers. It will be interesting to study by a mi- 
croscopic and totally self-consistent theory whether or 
not the center wave numbers of coupled antiferromag- 
netic spin fluctuations are close to (±1, ±1 ± 1/4) (tt/o) 
and (±1 ± 1/4, ±1) (tt/o) and whether those of coupled 
phonons are close to (0, ±1/2) (tt/o) and (±1/2, 0) (tt/o); 
those of antiferromagnetic spin fluctuations have been 
observed in neutron-scattering experiments and those of 
phonons correspond to the so-called 4x4 checkerboard. 
It has been proposed in a previous paper— that the co- 
existing state of antiferromagnetic moments and lattice 
distortion with such wave numbers is simply the observed 
4x4 checkerboard. As studied in AppendixO the coexis- 
tence of d-wave superconductivity and the antiferromag- 
netic state is possible, which implies that the coexistence 
of d-wave superconductivity and the checkerboard is also 
possible. 

The so-called zero-temperature pseudogap is observed 
in under-doped cuprate superconductors i^^i^^'^'^ A possi- 
ble scenario to explain it is the coexistence of da;2_j,2-wave 
superconductivity and another order parameter, which 
may be conventional or exotic. The large-gap phase stud- 
ied in 144. 41 is another possible scenario that can explain 
the zero-temperature pseudogap. The transition between 
the superconducting phase and the coexisting phase of 
superconductivity and the other order parameter must 
be a second-order transition because the symmetries of 
the two phases are different from each other, while the 
transition between the large-gap and small-gap phases is 
a first-order transition or a sharp crossover because the 
symmetries of the two phase are the same as each other. 
Nano-scale disorder or large inhomogeneity is observed in 
the gap structure of the pseudogap phase j^i^SiS which 
appears to support the scenario of a first-order transition 
or a sharp crossover rather than that of a second-order 
transition, i.e., the large-gap phase rather than the coex- 
isting phase. 

Another type of pseudogap is observed in the critical 
region above Tc- It is impossible to reproduce the pseu- 
dogap by the theory developed in this paper, which is 
restricted to T = K. In complete two dimensions, Tc 
decreases to -1-0 K because of critical fluctuationsi^ This 
fact implies that critical fluctuations are crucial in highly 
anisotropic quasi- two dimensions. It has been proposed 
in previous paper o^^'^^ that critical superconducting fluc- 
tuations are responsible for the pseudogap; the anisotropy 
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of the pseudogap is such that its size is proportional to 
r]i^{k) or [cos{kxa) — cos(fcya)]^, which is different from 
that of the superconducting gap. If the pseudogap opens, 
the single-site susceptibility Xs(0) is also enhanced even 
above Tc, as implied by results given in Appendix [Cl If 
the enhancement of Xs(0) is sufficiently large, antiferro- 
magnetic moments can appear; thus, the coexistence of 
antiferromagnetic moments and lattice distortion, which 
is simply the checkerboard, can also appear due to the 
electron-phonon interaction. The checkerboard has actu- 
ally been observed above T^j^l It will be interesting to 
examine whether or not and how the appearance of the 
checkerboard and the development of the pseudogap are 
correlated. Since the cuprate-oxide superconductor is in 
the vicinity of the Neel state, antiferromagnetic critical 
fluctuations are also so crucial that they are responsible 
for the so-called T-linear resistivity.— It is straightforward 
to include the effects of both superconducting and anti- 
ferromagnetic critical fluctuations within the framework 
of this paper. 



6. Conclusion 

In the vicinity of the Mott transition, a strong electron- 
phonon interaction arises from the modulation of the su- 
perexchange interaction by phonons. The Kondo-lattice 
theory of strong-coupling superconductivity based on the 
t-J model with the electron-phonon interaction is formu- 
lated. The self-energy of electrons is decomposed into 
the single-site and multisite self-energies. The calcula- 
tion of the single-site self-energy, which should be self- 
consistently determined with other properties such as 
the superconducting order parameter, the multisite self- 
energy, the phonon self-energy, the total exchange inter- 
action, and so forth, is reduced or mapped to a problem 
of self-consistently determining and solving the Anderson 
model. It is proved that the ground state of the mapped 
Anderson model is a conventional Fermi liquid, i.e., the 
single-site self-energy is that of a conventional Fermi liq- 
uid even if the order parameter is nonzero or the multisite 
self-energy is anomalous. The Fermi liquid characterized 
by the self-energy for the mapped Anderson model is fur- 
ther stabilized by the RVB mechanism. The density of 
states for electrons near the chemical potential is reduced 
by the RVB stabilization mechanism; it can be vanish- 
ingly small in a certain limiting case, although the den- 
sity of states for quasi-particles, which is directly related 
to the specific heat coefficient, is still nonzero and can 
be large. The stabilized Fermi liquid is a relevant unper- 
turbed state that can be used to study superconductivity 
and anomalous Fermi-liquid behaviors in the vicinity of 
the Mott transition or in the cuprate oxide. 

Not only the superexchange interaction, which arises 
from the virtual exchange of a pair excitations of electrons 
between the upper and lower Hubbard bands, but also the 
exchange interactions arising from that of a pair excita- 
tion of quasi-particles and that of a coupled excitation of 



spin fluctuations and phonons play a crucial role in the 
binding of dx^-yi-wetve Cooper pairs. On the basis of the 
analysis of the dynamical spin susceptibility, it is assumed 
in this paper that the imaginary part of the total ex- 
change interaction has a sharp peak or dip at ±uj* , where 
oj* ~ Wph in the normal state and iec ^ ^* ^ h^G +^ph 
in the superconducting state, where Wph is the energy 
of relevant phonons and ec is the superconducting gap. 
Then, it was shown that the dispersion relation of quasi- 
particles has kink structures at approximately ztw* above 
and below the chemical potential, the density of states for 
quasi-particles has dip-and-hump structures at approxi- 
mately ztw* outside the coherence peaks in the supercon- 
ducting state, and the anisotropy of the superconducting 
gap deviates from the simple anisotropy of (i3,2_j,2-wave 
superconductivity. These strong-coupling phenomena are 
consistent with observations in the cuprate-oxide super- 
conductor. As a functional of the exchange interaction 
or a function of the electron density, a first-order transi- 
tion or a sharp crossover can occur between small-gap and 
large-gap phases in the vicinity of antiferromagnetic insta- 
bility. It will be interesting to examine whether or not the 
large-gap phase is responsible for the zero-temperature 
pseudo-gap. 

In this paper, numerous phenomenological parameters 
and simplifications were introduced to simplify the nu- 
merical processes. It is desirable to develop a microscopic 
and totally self-consistent theory based on the formula- 
tion in this paper, in particular, to determine whether or 
not the imaginary parts of the spin susceptibility and the 
total exchange interaction actually have sharp peaks or 
dips at approximately ±u;*, as assumed in this paper. 



APPENDIX A: Proof of the Inequality of eq. ([3713]) 

In this Appendix, it is only assumed that E(j(e±iO, k)'s 
are analytical in the upper and lower half planes, respec- 
tively; they may or may not be singular. The inverse of 
the diagonal Green function is given by 

G-\e + iO, fe) = e + ^ - E{k) - E^(e + iO, k) 



|A,(e + iO,fc)|^ 



e- ^i + E{~k) + S_^(-e - iO. -fc) ' 

(Al) 



We define the following real functions: 

Si{e,k) = ReG^i(e + iO,fc) 

52(e,fc) =ImG;i(e + iO,fe) 



S^{e,k) 



Stick) + Siie,k)' 



^2"(e,fc) 



St{e,k) + Siie,ky 



(A2) 
(A3) 

(A4) 
(A5) 
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and 



52(e) = -Ini[r(e + iO) + S^(e + iO)] . (A6) 



In general, 



S2ie,k)>S2ie)>0, 



(A7) 



for any k. The single-site Green function is given by 
Giia{e + iO) — Yi{e) ~ iZi{e). According to the mapping 
condition p. lip . 

A(e) = -Imr(6 + iO) + ,,2,fl'l2,,^ ^ (^8) 



where 

X{e) = Z,ie) - S2ie)[Y,\e) + Zl{e)\. (A9) 
It is trivial that lo(e) = Zo{e), 

Y2{e) + Z^ie) = 1, (AlO) 

and 

Zi{e) > S2{e)Yo{e) = S2{e)Zo{e). (All) 
According to cqs. (|A10|) and (IA11|1 . it follows that 

Xie) = Zi(e) [Y2ie) + ^^(e)] - S2ie)[Y,\e) + Z^e)] 
>S2{e) [^Y,^{e)+Yo{e)Y2{e)] 

+ S2{e) [-Zfie) + Zo{e)Z2ie)] . (A12) 

Since the inequalities of 



Si{e,k) + Siie,k) 



and 



ly i- + S2ie,k)]' ^0 (^,4) 



7V^5f(e,fc) + 52^(e,fc) 



i.e., Yoie)x^ + 2Yi{e)x + Y2ie) > and Zo(e)a;2 + 2Zi(e)a; + 
Z2{e) > 0, hold for any real x, it follows that 



and 



Y,^{e)-Yo{e)Y2ie)<0, 



Z^e)- Zo{e)Z2ie) <0. 



(A15) 



(A16) 



According to eqs. (EI2), ([XTsjl . and ((XT6|l, i t 

follows that X{e) > 0. Thus, the inequahty ((3?T3)) . 
^A(e) > — Imr(e + iO), holds as a result of eq. (|A8|) even 
if the total self-energy So-(e ± iO, fc) is divergent or the 
order parameter A^{e ± iO, k) is nonzero. 



APPENDIX B: Dynamical Polarization Function in 
Spin Channels 

It is assumed that the dispersion relation of quasi- 
particles is given by 

^(fc) = —2tl [cos(fcj;a) + cos(fcj^a)] 

— 4t2COs{kxa) cos{kya) — n* , (Bl) 

where > 0, ^2 = —0.3^1, and /i* is such that 

(2/A^)Efe^[-?(fc)] = 0.85, i.e., ^i* = -1.0010|iJ|. Ac- 
cording to the Fermi-surface sum rule, the electron num- 
ber is 0.85 per unit cell in the normal state; it is also 
about 0.85 in the superconducting state. The Fermi sur- 
face or line for this set of parameters is similar to the 
observed one, as shown in Fig.[5Ka). The Green function 
for quasi-particles is given by 



e,(e±iO,fe) 



- m ± i7Q 
i77id(fe)AG 



-i?7id(fe)AG 
-I- e(fc) ± i7Q 



(B2) 



where rjid{k) = cosikxO) — cos(fcj,a). The energy depen- 
dence of the gap function Ag is ignored and the phase 
of Ag is chosen in such a way that Ag is positive. A 
small phenomenological lifetime width 7q = 0.01 |t|| is 
assumed for convenience in numerical processes if neces- 
sary. Figure [IKb) shows the density of states p*(e), which 
is defined by eq. (|4.25p . in the two cases of AG/|ti| = 
and Ag/I^iI — 1, i.e., in the normal and superconducting 
states; the superconducting gap is as large as 2Ag = 2|t5; | . 

Figure [TO] shows the polarization function P(a; + iO, q), 
which is defined by eq. p.36p . in the normal and super- 
conducting states for only three typical values of q near 
Qm = (±1, ±l)(7''/a) among many calculated values of 
q, which are distributed in the whole Brillouin zone. In 
the normal state, the real part of the static P(0, q) has 
a broad peak near Qm when it is plotted as a function 
of q. Since the peak is broad, the electron-phonon inter- 
action can play a crucial role in determining the center 
wave number Qsi of spin fluctuations or Qaf at which 
the static spin susceptibility x(0, q) is maximum. In the 
superconducting state, the imaginary part of P{u} + iO, q) 
is suppressed in the low-energy region of |w| < Aq by the 
opening of the gap and it has a peak near ~ Ag . The 
real part is enhanced in the low-energy region of w ~ Ag 
by the opening of the gap unless \q\a <C tt, which means 
that the opening of the superconducting gap assists the 
development of low-energy antiferromagnetic fluctuations 
and the appearance of antiferromagnetic moments. 

When the electron-phonon interaction is ignored, the 
spin susceptibility is described by 



1 



1^2 



[Po(t^ + iO,q)-P(c^ + iO,q)], 



(B3) 
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FIG. 9: (a) Fermi surfaces or lines in the normal state for the electron density (dashed line) n = O.li and (solid line) n = 
O.li + 0.05, where i is an integer, per unit cell. The thick solid line shows the Fermi surface for n = 0.85. (b) Density of states 
p*(e) for n = 0.85: (solid line) the normal state and (dashed line) the superconducting state. In the superconducting state, the 
logarithmic van Hove singularity is absorbed into the coherence peaks. 



(a) (iii) (b) (iii) 




FIG. 10: Polarization function P(a; + iO, q): (a) the normal state and (b) the superconducting state with IAg/^iI = 1; (i) 
q/{-K/a) = (0.8,0.8), (ii) (1,0.8), and (iii) (1, 1). Solid and dotted lines show the real and imaginary parts, respectively. 



where 

Po(c^ + iO,q) 



Xs(w + iO) 



Js{q) 



(B4) 



Following a previous study2S on the w-linear imaginary 
part of the susceptibility of the Anderson model, it can 
be shown that 



lim Im— 



1 



W2 x.(t^ + iO) 



P{lu + iO) 



0. (B5) 



Since the w-linear imaginary part vanishes, the energy 
dependence of Pq{uj + iO, q) is ignored in the following 
part. When a phenomenological parameter Hq is defined 
by Po(0, q) = P(0, q) [1 + Kg], it follows that 



X^(t^ + iO,g) 
X.(0,q) 



1 + K„ 



^(^ + io,g) 

^^(0,9) 



(B6) 



Figure [TT] shows + iO, <7)/xs(0, q) for q = Qm and 
four values of Kg in the normal and superconducting 



states. Because of the suppression of the imaginary part 
of P{uj + iO, q) in the superconducting state, the imagi- 
nary part of x(ti' + i0.q) shows a sharp peak near ut ~ Aq 
and a sharp dip near lu ~ — Aq. 

The susceptibility and the Green function for phonons 
are described in the spectral representation: 



p+oo 

Xsi'^^i,q) = / deAsq{e) 
Jo 



1 



ILUl 



iLUe 



+ 00 



1 



iuje — e iLUi + e 



and 

Jo 

At T = K, it follows that 

ksT'^Gxiii^e + iw„,q)xs(iwn,p) 

1 1 



1 



deAq^p{e) 



where 



Aq.pie) = / de'Axq{e~e')A,p{e'). 
Jo 



(B7) 



(B8) 



(B9) 



(BIO) 
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(a) (iv) (b) (iv) 




FIG. 11; Xs('^ + iOi 9)/Xs(0, q) for q — (±1, ±l)(7r/a): (a) the normal state and (b) the superconducting state with IAg/^iI ~ 1; 
(i) Kg — 0.8, (ii) 0.4, (iii) — 0.2, and (iv) 0.1. Solid and dashed lines show the real and imaginary parts, respectively. 



Note that eq. (|B9|) appears in Jsi.phi^^i,q) defined by 
eq. (|3.38p . As is shown in Fig. fTTT b). Asq{uj) with q — 
has sharp peaks near u = Aq in the superconducting 
state. Provided that Axq^Lo) has a sharp peak near uj = 
^ph, ^<j,p(e) also has a sharp peak near uj = Aq +Wph, so 
that Jsf-phli-^ + iO, q) has a sharp peak near lu — Ac+i-i^ph- 

In the complete theory, both Xs{^^i,q) and G\{iuji,q) 
should be self-consistently calculated with many other 
quantities. However, it is reasonable to expect that, in 
the final self-consistent solution, Imxs(w + iO,q ~ Qm) 
will have a sharp peak at oj ~ in the normal state 
and at w such that Aq ^ < Aq + Wph in the supercon- 
ducting state. 



superconductivity. 



APPENDIX C: Coexistence of Superconductivity 
and Antiferromagnetism 



According to eq. (|3.12p . it follows that 



P(e) = 



AA(e) 



{Rc[l/G.(6 + iO)]}%Ai(e) 



(CI) 



In general, Re[l/Gcr(-|-iO)] ^ for the non-half filling of 
electrons (n ^ 1). When a gap opens in p(e), a gap also 
opens in Aa(£). In this Appendix, the enhancement of 
the single-site spin susceptibility Xs(0) by the opening of 
the gap is studied. 

Although AA(e) should be self-consistently determined 
to satisfy the mapping condition ()3.11|) . the gap struc- 
ture of AA(e) is phenomenologically treated here. Since 
p{e) oc AA(e) oc |e| for small |e| at T = K for d-wave 



AA(e) 



An 



Ai, 

(Ai - Ao) 



Ao + (A2 - Ao) 

A2, 
0, 



-D <e< ~eo 
-eo < e < 

< e < eo 
\e\>D 



(C2) 

where eg ~ Aq and D ~ W/2, are assumed. According 
to eq. ((3lT|l . Ai ~ 0{W/tt) and Ai > A2 because the 
chemical potential is at the top of the lower Hubbard 
band. Although Ao = at T = K, Aq is treated as 
another parameter. It is assumed that Ai > Ao > 0. 

Since Uoo/D = +00, a doubly occupied configuration 
of d electrons is not allowed in any eigenstate. Thus, the 
lowest singlet state of the Anderson model is expanded in 
such a way that 

h(T 

+ Aka-p-.ciaCp-a + • • • ] |0) , (C3) 

where |0) is the Fermi vacuum for conduction electrons 
with no d electron. The ground-state wave function <I>s 
satisfies 



(7^A-AAAA+Hoxt)$s 



(C4) 



where Ha is defined by eq. (|2.22p . A/a — '^^nda + 

Y.ka Aa^ka, ^oxt = - Eo- (a/^ + CT^) n^a are iufinitcs- 
imally small external fields, and Eg is the energy of the 
singlet. When only Aq, Ad-ka, and Ak^.p^a are consid- 
ered, it follows that 



Es = / deAxie) 



Ed- AfL-e- Eg 
{Ed-Afi-e-Esf- 



(C5) 
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FIG. 12: l/7rAiXs(0) as a function of n: (i) Ao/Ai = 0, (ii) 
Ao/Ai = 0.1, and (iii) Ao/Ai = 0.2. Solid, dotted, dashed, 
dot-dashed, and double-dot-dashed lines represent eo/^i ~ 0, 
0.1, 0.2, 0.3, and 0.4, respectively; every solid line shows the 
relation of eq. (|C9|I . 



where Ed is the energy of the lowest doublet; 

A2 , D 

Ed = Ea + ed - n 1^ — , 

TT Al 



(C6) 



where Eq is the energy of the Fermi vacuum. 

It is easy to confirm that eq. (jCSp gives Eg such that 
Eg < Ed- When p{e) cx |e| for small |e|, the ground 
state of the Anderson model is a singlet; thus, the ground 
state of the t-J model within the restricted Hilbert sub- 
space where no order parameter exists is also a singlet, 
even if the reservoir effect is ignored. The d electron 
density, which is simply denoted by n here, is given by 
n = —dEs/dAfi, so that 



de- — 

D (Ed 



AA(e) 



EsY 



(C7) 



susceptibility is given by Xs(0) = [m/h^-^^^ or 



n (Ed-e- Es)^ 



(C8) 



When no gap structure develops in Aa(£), i.e., Ao = 
Al or £0 = 0, it follows that E^ = Ed - {2Ai/Tr){l-n)/n 
and 



X.(0) 



2Ai 1 



(C9) 



Since D 3> \Ed — Es\, D —> +00 is assumed here. When 
the gap develops in Aa(£), Xs(0) is enhanced or l/xs(0) 
is reduced, as shown in Fig. [121 

The Neel temperature of the t-J model is defined 

by 



-Js{q)--JQ{Q,q)--Jsi-plA^.q) 



The magnetization is given by m = —dEs/dh and the 



Xs{0) 4"-- 4" — ' 4 — -'-]^^^^ °' 

(CIO) 

where Tn should be maximized as a function of g. As 
shown in Fig. [T^l 1 /xs (0) is greatly reduced when a large 
gap opens, for example, when Aq/Ai ~ 0, Eq/Ai > 0.2, 
and n > 0.9. On the other hand, the superexchange in- 
teraction Js(g) is not reduced by the opening of the gap. 
As shown in Fig.fTUl P{0, q) is not reduced at least unless 
\q\a <C TT. According to eq. p.35p . P(0,q) is the main 
term of Jq(0, q). Thus, Jq (0,<7) is not reduced unless 
\q\a ^ TT, which implies that Jsf-ph(0, g) is also not re- 
duced unless |<7|a ^ tt. If a superconducting gap 2Ag as 
large as 2 Aq/Ai ^ 0.4 opens, it is probable that antifer- 
romagnetism appears for n > 0.9, at least at T = K, 
i.e., antiferromagnetism coexists with superconductivity 
for n > 0.9. 
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